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Abstract. We consider a relativistic quantum model of confined 
'. massive spinning quarks and antiquarks which describes leading 

Regge trajectories of mesons. The quarks are described by the 
J> ; Dirac equations and the gluon contribution is approximated by 

the Nambu-Goto straight-line string. The string tension and the 
current quark masses are the main parameters of the model. Addi- 
tional parameters are phenomenological constants which approxi- 
qq \ mate nonstring short-range contributions. Comparison of the mea- 

sured meson masses with the model predictions allows one to de- 
^ ' termine the current quark masses (in MeV) to be m s = 227±5, m c = 

1440 ± 10, m b = 4715 ± 20. The chiral SU 3 model[23] makes it 
Ph! possible to estimate from here the u- and rf-quark masses to be 

jjn ; m M = 6.2±0.2 Mev and m d = 11.1 ± 0.4 Mev. 



PACS Numbers: 12.15. F, 12.50.C, 11.17 



It has been believed for a long time that properties of quarks 
confined in a meson are closely related to those of the relativis- 
tic string with Nambu-Goto self-interaction [1]. The anomaly in 
the quantum string theory in 4-dimensional space-time has led to 
other important applications of the string theory [1]. Nevertheless, 
the hadron theory can use particular simple configurations of the 
string for approximate description of the hadrons if these configu- 
rations admit relativistic quantization. If the approximate hadron 
model obtained in this way appears to be in acceptable agreement 
with the experiment one can try next more complicated string 
configuration, having in mind that at some step the whole notion 
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of string may fail, especially when more experimental information 
about hadron daughter trajectories will be available. 

The simplest string configuration, a straight-line string, was 
quantized in [2,3] in accordance with the Poincare invariance and 
gave good agreement with the spectrum of the light-quark mesons 
lying on the leading Regge trajectory. The next approximation 
was to take into account the masses and the spins of the quarks 
attached to the ends of the string. This has been done in [4-16] 
with various assumptions. 

The distinctive features of the present approach as compared 
with those of refs[4-16] are the consistent treatment of the quark 
spins and the canonical quantization. The gauge invariant formal- 
ism is used throughout the paper. We also show that there is no 
radial motion of the quarks along the rotating straight-line string. 
This means that the daughter meson states correspond to higher 
modes of the string (vibrations). 

The advantage of the present approach as compared with the 
potential models ([17] for example) is relativistic invariance (in 
[17] it is only approximate) and use of current quark masses (in 
[17] constituent quark masses are used). The disadvantage of the 
present paper is restriction to the leading Regge trajectories, i.e., in 
the potential model language, to the lowest radial excitations. 

We consider the Nambu-Goto straight-line string with point-like 
massive spinning quarks attached to its ends. This is an extended 
relativistic object [18,19] called rotator for which the explicitly rel- 
ativistic description introduces auxiliary variables resulting in a 
symmetry of the rotator Lagrangian. The Hamiltonian of the ro- 
tator is given by an implicit function which can be calculated nu- 
merically. For important particular cases ( light or heavy quarks) 
series expansions for the Hamiltonian are obtained. 

The quark spins are described by anticommuting spin variables 
obeying constraints [20]. Special care has been taken to ensure 
conservation of these constraints [21,22]. 

Canonical quantization of this system preserving the Poincare 
invariance yields meson states with different spins and parities lying 
on Regge trajectories which depend on the quark masses. The 16- 
component wave function of a composite meson satisfies two Dirac 



2 



equations and a spectral condition which can be compared with 
the experimental mass spectrum. 

The spectral condition contains a contribution of the universal 
string confining mechanism together with nonstring short-range 
contribution which is treated phenomenologically. The dominant 
part of the short-range contribution do not depend on the meson 
spin J and its decreasing with J part is seen only in low- J quarko- 
nia. The string contribution dominates when at least one quark is 
light and grows with the meson spin. On the other hand, it is near 
threshold for low-spin heavy-quark mesons. The string contribu- 
tion to the T(1S) mass is about 200 MeV and to the x 62 (l.P) mass 
- 350 MeV. 

So, the present approach in its simple form is applicable to 
mesons containing at least one light quark where the nonrelativis- 
tic potential models are not applicable. For heavy quarkonia the 
string mechanism should be supplemented with other small ( com- 
pared to heavy-quark masses ) contributions to account for the fine 
structure of the levels. 

We compare the model with experiment for the trajectories 
with P = C = (— 1) J and lowest states having J PC = 1 . For these 
trajectories mesons with highest spins were observed and mixing 
with other trajectories is negligible. 

This comparison with meson masses allows to estimate the cur- 
rent s-, c- and 6-quark masses assuming that the current u- and 
d-quark masses are zero within error bars. We then use the chiral 
SUs model [23] to estimate the u- and c?-quark masses through the 
s- quark mass to check the consistency of the calculations. 

To check the model we have used the obtained quark masses to 
calculate the masses of mesons not used in the input. We compare 
the predicted masses with experiment and with the results of the 
potential model [17] and discuss a possible interpretation of the 
gluon string in terms of the potential model. 

So, let us consider a simplest extended relativistic object - a 
straight-line 

x(r,a) = r(r) + f(r,a)q(r), (1) 
where r is a 4-vector corresponding to a point on the straight-line, 
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q is an affine 4-vector of its direction, / is a scalar monotonic func- 
tion of <t labelling points on the line and r is a scalar evolution 
parameter. We shall not fix the coordinates /j(r) = /(r, <7j(r)) of the 
end points of the string considering them as dynamical variables 
to be determined from extremum of an action. Then the explicit 
Poincare covariance of (1) introduces superfluous variables not nec- 
essary for description of the straight-line as a physical object, so 
that theory in terms of (1) must be invariant under a group of 
three sets of r-dependent transformations ( gauge transformations 

) 

1) shift of r along q: 

r^r + f(r)q, (2) 

2) multiplication of q by an arbitrary scalar function: 

q^g(r)q, (3) 

3) reparametrization of r, what means that the Lagrangian must 
satisfy the condition 

C(h(T)z,h(r)(h(r)z)- ) = h(r)C(z,z), (4) 

where z and z mean every r-derivative in the Lagrangian. 

This symmetry implies that the phase-space variables of our 
system obey three constraints which are in involution with respect 
to their Poisson brackets; the canonical Hamiltonian is zero and 
the total Hamiltonian is a linear combination of the constraint 
functions. 

Invariants of a symmetry play an important role in description 
of a symmetric system. In our case they are orthonormal vectors 
along line direction, velocity of the line rotation and velocity of its 
movement as a whole 

n=(-q 2 )' 1/2 q, v 1 =b~ 1 h, v° = (rl)~ 1/2 r ± , (5) 

where 

b = {-n 2 ) 1 ' 2 (6) 
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and 

f- k ± = (g kl + n k n l + v lh v ll )fi. (7) 

The angular velocity b is invariant under (2) and (3) and transforms 
as the Lagrangian under (4). The scalar invariant of the symmetry 
is 

/ = r 1 (fi) 1 / 2 . (8) 

We shall label points on the string with respect to the instant 
center of its rotation z 

(-q 2 ) 1/2 f = z + y, (9) 

z = b'Hv 1 (10) 

(velocity of the point r + zn, orthogonal to q, is orthogonal to v 1 
). The length of the rotator at fixed r is \y 2 — Vi\- From x\ > it 
follows that \yi\ < I. 

We shall take quark spins into account later on. Without quark 
spins the Lagrangian of our model is a sum of the Nambu-Goto 
Lagrangian for an open string with a string-tension parameter a 
and two Lagrangians for free point-like particles with masses mi 
and m 2 and velocities of the ends of the string 

C = -a I" g^da -Y,™^) 1 ' 2 , (11) 

J<T1 i 

where g = (xx') 2 —x 2 x' 2 is minus determinant of the induced metric of 
the string worldsheet and ii = dx(r, <7j(r)) / dr, i = 1,2 are velocities 
of the string ends. Using the notations introduced above we can 
rewrite (11) for the straight-line string (1,9) in the form 

C = -bF, (12) 

where F is a gauge and Poincare invariant function 

F = a [ y \f - x 2 f/ 2 dx + Y,mS 2 ~ vt ~ (13) 

w i = b~ 1 (y i + z -rri). (14) 
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We shall consider the case when 

by^O (15) 

(this is a gauge invariant condition). Then we must consider Wi as 
independent variables and the stationary condition with respect to 
them yields 

Wi = 0. (16) 

The other way to obtain this result [16] is to consider the Euler- 
Lagrange equations following directly from (11) which give for the 
straight-line string 

y — 0, z — rn = 0. (17) 

Eq.(16) follows from here by continuity. 
We conclude that for our model 

F = a( y \l 2 -x 2 ) l / 2 dx + Y.m t {l 2 -y*y /2 (18) 

with yi satisfying the stationary condition 

dF/dyi = 0, (19) 

or 

(-iy yi = (l 2 + K/2a) 2 ) 1/2 - (m t /2a). (20) 
Calculating the momenta p and tt canonically conjugate to r and 

Q 

p = -dC/df, tt = -dC/dq (21) 

we get three constraints 0j = 0, i = 1,2,3 where the constraint 
functions are 

0i = pq, 02 = Trq, (22) 
03 = L - K. (23) 

Here 

L = ((g 2 - M7pV) 1/2 (24) 

is the magnitude of the conserved orbital spin 

Lfj, = e^ pa p"M^/2m (25) 
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where 

AP" = r V 1 + q [fl K u] (26) 

is the angular momentum tensor. K is a function of m = (p 2 ) 1 ^ 2 , 
implicitly given by the equations 

K = lm-F, (27) 

dF/dl = m. (28) 

The rotator Hamiltonian is a linear combination of the con- 
straint functions 

H= Yl ( 29 ) 

i=l,2,3 

It determines the dynamical equations for any variable X 

X = {X,H}, (30) 

(pi = after calculating the brackets and the non-zero Poisson brack- 
ets are 

{jry} = {<fX} = gr. (31) 
We can choose gauge conditions to fix c\ $ = in (29): 

pre = 0, q 2 + 1 = 0. (32) 

To obtain the Poisson brackets in this gauge we introduce new 
variables having vanishing brackets with the constraints (22) and 
(32) 

p, r = r + {{prr)q - {pq)n)/p 2 , v = (-q 2 p ) 1/2 q p , L (33) 

(ip = {d^ u ~~ V il P v /p 2 )Qv)- To have zero brackets of the external coor- 
dinate of the rotation center r with the internal coordinates v and 
L we use four orthonormal vectors e a , a — 0, 1, 2, 3 

e = p/m, e a e/3 = g a p (34) 

and introduce new variables 

n a = -e a v, L a = -e a L, a =1,2,3, (35) 
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1 de u 

z = r + -e a b c e a u-r^L c . (36) 
The non-zero Poisson brackets of the new variables are 

{p k ,z l } = g kl , {L a ,L b } = e abc L c , {L a , n b } = e abc n c . (37) 
The constraint function 3 now takes the form 

3 =((£ a ) 2 ) 1/2 -^(m) (38) 

and the solution of the dynamical equations (30) can be easily 
obtained to be 

z = z + IVp/m, (39) 

n = n cos V — rix sin V, (40) 



V = J c 3 dr. (41) 
Prom (39) the laboratory time of the rotation center 

t = z° - z° = lVp°/m (42) 
and the space coordinates of this point 

z a = z% + p a t/p° (43) 

correspond to its movement in the laboratory with constant veloc- 
ity p a /p°. The direction of the rotator rotates with constant angular 
velocity 

- - < 44 > 

where / = l(m) from (28). 

The canonical quantization can now be performed quite easily. 
We replace our variables by operators and their Poisson brackets 
(37) by commutators. The constraint equation now holds for the 
wave function 

[((L a f) 1 / 2 -K(m)-a }^ = 0, (45) 

where in the operator form of (38) we have added a term ao to 
account for nonstring short-range contributions. 
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Our quantum system is relativistic because the quantization 
procedure transforms the classical Poisson brackets of and M ua 
into commutators without any change in their form, so that the 
Poincare algebra is fully preserved. 

Quark spins are important especially for small L. They were 
taken into account in [21,22] where the spinless-particle Lagrangians 
in eq. (11) were replaced by those of Berezin and Marinov [20] and 
a special term was added to preserve conservation of the spin con- 
straints, with the result that for the leading Regge trajectories one 
can simply replace the orbital spin L in (45) by the total meson 
spin J. This yields 



for the physical eigenstates with fixed dependence of space and 
charge-conjugation parities P and C on J. 

The function K{m) is given by eqs. (18,20,27,28). We must solve 
eq.(28) to find I as a function of m and put this function into (27). 
This can be done numerically for any quark masses. For important 
particular cases K can be expanded into series. For light quarks 



(J(J+1)) 



1/2 



K{m) + a 



(46) 





For heavy quarks 

D = m — mi — m 2 <C Wj (49) 



A . (m) .i ( |^/v' 2 (i4^ + o ( (^n), (so. 



For light and heavy quarks 

, ixmi Id 
d = m-m 2 , yi = -— < 1, x 2 = < 1, (52) 



nm 2 



K(m) 




2tt 2 



35 




(53) 
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We see that the slope of the trajectory for mesons formed by a 
heavy and a light quark (antiquark) is twice as big as for light- 
quark mesons. 

The term a in (46) can in general depend on J, but it can not 
grow with J. An analysis of Coulomb-like short-range interaction 
suggests the following dependence of a on J ( or on m i what is 
practically the same when eq.(46) is fulfilled) 



where A and B do not depend on J. In all cases considered be- 
low the first term in (54) dominates, so the precise form of the 
second term is not important for our conclusions. As a first ap- 
proximation one could neglect the second term and to get the quark 
masses within error bars following from comparison with experi- 
ment. On the other hand the second term allows one to get good 
agreement with the experimental heavy-quarkonia spectrum. The 
errors in the quark masses in this case formally reduce and to es- 
timate their values one have to go outside of the model and to 
analyse the interaction between mesons and their decay channels. 
An approximate analysis of this problem was performed in ref. [17] 
with the result that the error in the heavy- quark meson masses 
is about 10 Mev. We tentatively take this value as an error in 
the heavy-quark masses deduced from a precise fit to experimental 
meson masses with the help of the second term in (54). 

Assuming B in (54) to be of order 1 we see that the second 
term in (54) is negligible when one or both quarks are light. It is 
negligible also for the ss-mesons below. 

We shall apply Eq.(46) to the leading trajectories with P = C = 
(— 1) J and the lowest states having J PC = 1 . Estimates show they 
do not mix with other trajectories with the same J PC having much 
heavier states. 

Applying eq.(46) to the leading p - and ^-trajectories we have 




(54) 



K{m pJ ) = K(m K *j), 



(55) 
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or, neglecting the u- and d-quark masses 



A/J\l+'r + -L Zj + (/f)), (56) 
rriK*.j tc 10 lo7r 

3vr m 2 j 
ZJ = T (1 " m^~ ] - (57) 

The error from neglecting the u- and d-quark masses can be esti- 
mated from the cu- and p-mass difference to be 1.8%. Using exper- 
imental data for the meson masses from [24] we obtain the corre- 
sponding values for the strange quark mass shown in Table 1. The 
error in the average m s corresponds to the accuracy of calculations 
and, partly, to the accuracy of the model. 

We get the following values for the other model parameters 

a = .176, 2na = a'' 1 = 1.11 GeV 2 , (58) 

a = A = .88. (59) 

The parameter (59) is the same for the light and the strange quarks 
and corresponds to the intercept parameter (of J with the K = 
axis) J = .51. 

Knowing the strange-quark mass we can estimate the light- 
quark masses from the linear approximation of the chiral SU3 model 
[23] 

m u /m d = 0.554 ±0.002, m s /m d = 20.13 ± 0.03. (60) 
Using here m s from Table 1 we get (in Mev) 

m M = 6.2±0.2, m d = 11.1 ± 0.4. (61) 

We see that neglecting these masses in the above calculations does 
not introduce any noticeable error. 

To check these results we can use them to calculate masses of 
mesons consisting of ss, Table 2. They are in good agreement with 
the experimental values. 

To obtain the c-quark mass we consider eqs. (46), (54) for the 
D* and D2 _mesons - The second term in (54) is negligible and 

y/2 = K(D*)+A(c), V6 = K (D*) + A(c) , (62) 
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what allows us to calculate the c-quark mass through those of D* 
and D\ (Table I) and to estimate A(c): 

A(c) = .90. (63) 

We see that it is close to the constant A for the light quarks (59). 
To describe this closeness let us remark that the shift 0.02 in a 
yields the shift from —10 to —17 Mev in the vector-meson masses. 
This shift decreases for higher meson spins. 

Application of eqs.(46),(54) to the cc-mesons J/ip and ^(lP) 
gives the constants 

A(cc) = 0.90, (64) 
what coincides with (63) , and 

B(cc) = 1.43. (65) 

For the 6-quark we can not carry out a similar analysis because 
the mass of B\ is not known. To get an estimate of the 6-quark 
mass we have to rely on an assumption. The safest assumption 
seems to be 

A(b) = A(bb) (66) 

similar to the case of the c-quark (63), (64). Using the masses of 
B*,T(1S) and ^(l^O-mesons we get the 6-quark mass in Table I 
and 

A(b) = A{bb) = 0.77, (67) 

B{bb) = 3.14 (68) 

Experimental measurement of the B\ mass is important for check- 
ing the assumption (66). 

Now we can calculate masses of other mesons belonging to 
our trajectory. Some of them are presented in Table 2, together 
with experimental data available and predictions of the potential 
model of ref.[17]. This model is based upon linear rising poten- 
tial, Coulomb-like short-range potential from perturbative QCD, 
approximate relativistic corrections and constituent quark masses 
among other parameters. 
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It is tempting to conclude from Table 2 that the present model 
slightly better agrees with the data and that future precise mea- 
surements might distinguish both models. But far more impressive 
is the similarity of the results of apparently quite different calcu- 
lations. This similarity confirms the main physical motivation for 
considering the gluon string, namely, the string describes two sepa- 
rate mechanisms of the potential approach, confining potential and 
the constituent quark masses. 

In conclusion, let us discuss the relation between quark masses 
in this model and in QCD. The present model is a quantum me- 
chanical model of free quarks bound in mesons. Since it agrees 
with experimental data it is reasonable to assume that the quark 
masses of this model are the current quark masses entering as pa- 
rameters into the QCD Lagrangian when one uses the on-mass-shell 
perturbative renormalization procedure summed to all orders. 

It would be interesting to check the obtained values of the cur- 
rent quark masses in other applications. 

The author is grateful to V. A.Petrov, Yu.F.Pirogov and A.V.Razumov 
for discussions and to Prof. A. D. Krisch for the kind hospitality 
at the University of Michigan where this work was finished. 

Table 1: The input meson masses and the predicted current quark masses in 
the present model. 



Meson 


Input meson 


Quark masses in MeV 


spin J 


masses [24] 


Calculated in 
the present model 


Other estimates [24] 


1 


p, K* 


m s = 220 ± 4 




2 




m s = 234 ± 4 




3 


Ps, Kt 


m s = 204 ± 18 






average 


m s = 227 ± 5 


m s =100 to 300 


1 


D*,D* 2 


m c = 1440 ± 10 


m c =1.0 to 1.6 GeV 


1 


T, B*, X b2 


m b = 4715 ± 20 


m b =4.1 to 4.5 GeV 



References 
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Table 2: The model predictions for meson masses (in Mev) and comparison 
with the potential model predictions of ref. [17] (q stands for u or d). 



Quark 


Meson 


Present 


Experimental 


Potential 


content 


spin «r ° 


model 


values 


model [17] 


qq 


2++ 


1317 


1318.1 ±0.7 


1310 




3" 


1690 


1691 ± 5 


1680 




4++ 


1993 




2010 




c 

5 


2255 




2300 


qs 


A + 


OAOA 

2UoU 


on a c _i_ n 


2111) 


ss 


i — 
1 




1 m n /i 1 o _i_ n nno 


1U2U 




o++ 
z 






lOOU 




3~ 


1873 


1854 ±7 


1900 




4 ++ 


2160 




2200 


cq 


3" 


2780 




2830 


cs 


\~ 


2134 


2112.4 ±0.7 


2130 




2+ 


2561 


2573.5 ± 1.7 


2590 




3~ 


2870 




2920 


cc 


3— " 


3830 




3850 


bq 


2+ 


5720 




5800 


bs 


\- 


5430 




5450 


bc 


\- 


6410 




6340 


bb 


3~ 


10110 




10160 
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